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Abstract. The complete intcgrability of a generalized Riemann type hydrodynamic system 
is studied by means of a novel combination of symplectic and differential-algebraic tools. A 
compatible pair of polynomial Poissonian structures, a Lax representation and a related infinite 
hierarchy of conservation laws are constructed. 



1. Introduction 



o 

-t— > 
o 

O 

r -j- ■ We shall study the complete integrability of the dispersionless Riemann type hydrodynamic flow 

(1.1) £>f = z 2 x , D t z = 

on a 27r-periodic functional manifold M N C C (oo) (R/2ttZ; R n ), where N G N is an arbitrary 
natural number, the vector (u, D t u, D 2 u, Df~ 1 u, z) T G M N , the differentiations D x := d/dx, 
D t := d/dt + ud/dx satisfy the Lie-algebraic commutator relationship 



(1.2) [D x ,D t ]=u x D x . 



and t G R is an evolution parameter. The system can be considered as a slight generalization of 
the dispersionless Riemann hydrodynamic system (suggested recently by M. Pavlov and D. Holm 
' [S]) in the form 

(1.3) D?~ x u = z, D t z = 

^ . 

for N G N and extensively studied in [lj [3l SI [5j [21 [6] , where it was proved that it is a Lax integrable 
bi-Hamiltonian flow on the manifold M N and possesses an infinite hierarchy of mutually commuting 
dispersive Lax integrable Hamiltonian flows. 

For the case N = 2 it is well known [8j [11] that the system (jl.ip is a smooth Lax integrable 
bi-Hamiltonian flow on the 27r-periodic functional manifold M 2 , whose Lax representation is given 
, by the compatible linear system 



(L4) DJ ~ ( -\(u + u x /z x ) -z xx /z x ) /; Dtf - ( -Xz x u x ) ) /! 

where / G C^(R 2 ; R 2 ) and AeR is an arbitrary spectral parameter. 

Our focus here is an investigation of the Lax integrability of the Riemann type hydrodynamic 
system JTIJ) for N = 3 on a 2?r-periodic functional manifold M 3 C C (oo) (R/2ttZ; R 3 ) for a vector 
(u, v,z) T G M 3 . We treat this problem in the following extended form: 

(1.5) D t u = v, D t v = z 2 x , D t z = 0. 

The flow (11.51) can be recast as a one on a 2n- periodic functional manifold M 3 C (R/2ttZ; 
for a vector (u, v, z) T G M 3 as 

(1.6) D t u = v, D t v = z, D t z = —2zu x , 
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where, for further convenience, we have made the change of variables: z :- 
the form of the flow (|1.6[) : 

(1.7) 



We will also use 




K[u, v, z] 



v — uu x 
z — uv x 



defining a standard smooth dynamical system on the infinite-dimensional functional manifold M 3 , 
where K : M 3 — > T(M 3 ) is the corresponding smooth vector field on M 3 . 

In the sequel, we shall prove the following result using symplectic gradient-holonomic and dif- 
ferential algebraic tools. 

Proposition 1.1. The Riemann type hydrodynamic flow |J.7| ) is a bi-Hamiltonian dynamical 
system on the functional manifold M 3 with respect to two compatible Poissonian structures d, rj : 
T*(M 3 ) T(M 3 ) 



(1.8) 




2z 1 / 2 D r z 1 / 2 



d- 1 Uxd- 1 

, ?y := ( d~ x u x Vxd- 1 + d~ 1 v x d~ 1 z x -2z 
z x d~ l +2z 



possessing an infinite hierarchy of mutually commuting conservation laws and a non- autonomous 
Lax representation of the form 



(1.9) 



Dtf = 



D x f = 




\Vyfz 

-XHv^/z 

-\v x /sfz+ 



A 3 (to 2 



z 

-Xtz 



A 2 v / i(u — tv)- 
-z x /2z 



f, 



where A € R is an arbitrary spectral parameter and f G C^^R 2 ;] 



2. Symplectic gradient-holonomic integrability analysis 

Our first steps in proving Proposition 1.1 are fashioned using the symplectic gradient-holonomic 
method, which takes us a long way towards the desired result. 

2.1. Poissonian structure analysis on the functional manifold M 3 . By employing the sym- 
plectic gradient-holonomic approach [51 [TTJ [T^] to studying the integrability of smooth nonlinear 
dynamical systems on functional manifolds, one can find a set of conservation laws for (|1.7[) by 
constructing some solutions tp := ip[u, v, z] G T*(M 3 ) to the functional Lax gradient equation: 

(2.1) dtp/dt + K'> V = grad£, 

where <// = tp'>*, C £ D(M 3 ) is a suitable Lagrangian functional and the linear operator K''* : 
T*(M 3 ) T*(M 3 ) is the adjoint with respect to the standard convolution (•,•) on T*(M 3 ) x 
T(M 3 ), of the Frechet-derivative of a nonlinear mapping K : M 3 —> T(M 3 ); namely, 

(uD x -v x 
1 u x + uD x 
1 -u x 

The Lax gradient equation (|2.ip can be, owing to (|1.3p . rewritten as 

(2.3) D t ip + k[u, v, z](p = grad £, 
where the matrix operator 

( z x 

(2.4) k[u,v,z] 
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The first vector elements 

(2.5) ip§ [it, v, z] — (z — uv x , —v + uu x , u),C$ =0 
<P v [u,v,z] = (v x ,-u x ,-l) J ,£ v = 0, 

if [u,v,z] = {-(u x z-V%,(z- 1 /%,(v x /2-u 2 x /4)z- 3 / 2 y,£ = 0, 

as can be easily checked, are solutions of the functional equation (|2.3p . From an application of 
the standard Volterra homotopy formula 

(2.6) H := / d(j,((p[[iu, fiv, fj,z], (u, v, z) T ), 

Jo 

one finds the conservation laws for (|1.3|) : namely, 



(2.7) H ri = - I dx(2uz - v A - u 2 v x ), 



o 

H# := / dx(uv x /2 - vu x /2- z), H := - / dx{u 2 x - 2v x )z~ 1/2 . 
Jo * Jo 

It is now quite easy, making use of the conservation laws (|2.7[) . to construct a Poissonian 
structure i? : T*(M 3 ) T(M 3 ) for the dynamical system (fT?) . If we use the representations 

(2.8) H# = / dx(uv x /2~ vu x /2- z) :^ (^,(u x ,v x ,z x y), 

Jo 

^ := (-u/2, U /2,z- 1/2 D-V/ 2 /2) T , 
it follows that the vector -0^ € T*(M 3 ) satisfies the Lax gradient equation (|2 . 3|) : 

(2.9) D t ip# + k[u, v, z]ip$ = grad 

where the Lagrangian function £,j = (^,-ftT) — H$. Thus, based on the inverse co-symplectic 
functional expression 

/ -1 

(2.10) d- 1 : =^-V^*= 1 

V z- 1 / 2 D x 1 z- 1 / 2 /2 

one readily obtains the linear co-symplectic operator on the manifold M 3 : 

1 

(2.11) i?:= | -1 

2z 1 / 2 D x z 1 / 2 

which is the corresponding Poissonian operator for the dynamical system (jl.2|) . It is also important 
to observe that the dynamical system (|1.2|) is a Hamiltonian flow on the functional manifold M 3 
with respect to the Poissonian structure (|2.11j) . 

(2.12) K[u, v, z] = grad H n . 

2.2. Poissonian structure analysis on M 3 . In what follows, we shall find it convenient to 
construct we will construct other Poissonian structures for dynamical system (11.5[) on the manifold 
M 3 , rewritten in the equivalent form 





(2.13) — \ v \ = K[u,v,z 



where K : M 3 — > T(M 3 ) is the corresponding vector field on M 3 . To proceed, we need to obtain 
additional solutions to the Lax gradient equation (|2 . 3[) on the functional manifold M 3 



(2.14) D t ip + k[u, v, z]ip = grad C, 
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where the matrix operator is 

(2.15) k[u,v,z]:= 

and which we may rewrite in the componentwise form 

(2.16) D t 4> {1) = v x 4> {2) + z x ^ 3) + SC/Su, 

D t 4> {2) = - u x 4> {2) + SC/Sv, 

D t i>W = 2(z x ^) x - u x ^ 3) + SC/Sz, 

where the vector -tp := { r 4> < ^'\ "4 ! '"' 2 \ ^ 3 ^) T £ T*(M 3 ). As a simple consequence of (|2.16[) . one obtains 
the following system of differential relationships: 

£>3^(2) = _ 2 f2^(2) +D*d- 1 (6C/5v)- 
/ 917 x -d^ 1 < grad C, (u x ,v x ,z x ) J >, 

[ D t jW = -^ + d- 1 (S£/Sv), 

Dtfl*) =2z x ^ x 2) +d- 1 (SC/6z). 

Here we have defined ^ (2) , ^ (3) ) T := , ^ , </>£ 3) ) T and made use of the commutator 

relationship for differentiations D t and D x : 

(2.18) [D t ,a- 1 D x ]=0, 

which holds for the function a := l/z x , where D t z = O.It therefore follows that after solving the 
first equation of system (|2.17[) . one can recursively sole the remaining two equations. In particular, 
it is easy to see that the three vector elements 

^ = (-v,u, -2z x ) J , £ = 0; 

(2.19) j> B = {-u x /z x , l/z x , (u 2 x - 2v x )/(2z 2 x ))\C e = 0; 

^ = {u/2,-x/2,d- 1 [(2v x -ul)/{2z x )]),C n = (D^K) - H#, 

are solutions of the system ()2.17j) . The first two elements of (|2.19|) lead to the Volterra symmetric 
vectors tpo = D x tp ,ipg = D x ?pg £ T*{M 3 ) : -tp' = $q* , = entailing the trivial conservation 
laws (ipo,K) = = (tp0,K). The third element of (|2.19|) gives rise to the Volterra asymmetric 
vector ip v := D x ipr) '■ ipL 7^ V'n*) entailing the following inverse co-symplectic functional expression: 

/ d -df- \ 

df 



(2.20) fj- 1 := - $f = 



25 



Correspondingly, the Poissonian operator r\ : T*(M 3 ) — > T(M 3 ) is 

/ d- 1 

(2.21) fj = O- 1 ^ v x d- 



d- 1 u^- 1 



z x d- 1 

subject to which the following Hamiltonian representation 

(2.22) K[u,v,z] = -fj grad H^U= Z 2 
holds on the manifold M 3 . 

2.3. Hamiltonian integrability analysis. Next, we return to our integrability analysis of the 
dynamical system (|1.7|) on the functional manifold M 3 . It is easy to recalculate the form of the 
Poissonian operator (|2.21|) on the manifold M 3 to that acting on the manifold M 3 , giving rise to 
the second Hamiltonian representation of (|1.7[) : 

(2.23) K[u, v, z] = -7] grad H#, 

where rj : T*(M 3 ) — > T(M 3 ) is the corresponding Poissonian operator. As a next important point, 
the Poissonian operators (|2.11j) and (|2.21j) are compatible [TDl El El [7] on the manifold M 3 ; that 
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is, the operator pencil (i9 + Xtj) : T*(M 3 ) — > T(M 3 ) is also Poissonian for arbitrary A G R. As a 
consequence, any operator of the form 

(2.24) & n := ^iT 1 ??)" 

for all n G Z is Poissonian on the manifold Af 3 . Using now the homotopy formula (|2.6p and 
recursion property of the Poissonian pair (|2.12j) and (|2.21j) . it is easy to construct the related 
infinite hierarchy of mutually commuting conservation laws 

( 2 25) = Jo dfi(gTad jj[(J,u, fiv, fiz], (u, v, z) T ), 

grad 7j [w, v, z] :— A J grad H v , 

for the dynamical system (|1.7[) . where j G Z + and A := d^ lr q : T*(M 3 ) — > T*(M 3 ) is the corre- 
sponding recursion operator, which satisfies the so called associated Lax commutator relationship 

(2.26) dk/dt = [A, K' 1 *]. 

In the course of above analysis and observations, we have proved the following result. 

Proposition 2.1. The Riemann hydrodynamic system \1. ?] ) is a bi-Hamiltonian dynamical system 
on the functional manifold M 3 with respect to the compatible Poissonian structures i),r/ : T*(M 3 ) 
T(M 3 ) 

d- 1 Uxd" 1 




(2.27) t? := I —1 1 77 ;= I d^u* v^' 1 + d- x v x d- 1 z x -2z 

z 1 ' 2 j \ z x d~ 1 +2z 

and possesses an infinite hierarchy of mutually commuting conservation laws (2.25\) . 

Concerning the existence of an additional infinite and parametrically R 3 A-ordered hierarchy 
of conservation laws for the dynamical system (|1.2j) . it is instructive to consider the dispersive 
nonlinear dynamical system 

/ du/dr \ / -(z- 1/2 )x \ 

(2.28) dv/dr = -i? gr&dH [u,v,z] := -{u^- 1 / 2 ),, \=K[u,v,z\. 
V dz/dr J \ z V*{±g±) x J 

By solving the corresponding Lax equation 

(2.29) d<p/dt + K''*<f = 

for an element (p £ T*(M 3 ) in a suitably chosen asymptotic form, one can construct an infinite 
ordered hierarchy of conservation laws for (|1.2[) . which we will not delve into here. This hierarchy 
and the existence of an infinite and parametrically R 3 A-ordered hierarchy of conservation laws for 
the Riemann type dynamical system (|1.2|) provided compelling indications that it is completely 
integrable in the sense of Lax on the functional manifold M 3 . We will be the complete integrability 
in the next section using rather powerful differential-algebraic tools that were devised recently in 

una©. 

3. Differential- algebraic integrability analysis: N = 3 

Consider a polynomial differential ring JC{u] c JC := R{{a;,£}} generated by a fixed functional 
variable u £ R{{x,t}} and invariant with respect to two differentiations D x := d/dx and D t := 
d/dt + ud/dx that satisfy the Lie- algebraic commutator relationship (jl.2p 



(3.1) [D x ,D t ]=u x D x 

together with the constraint (|1.6|) expressed in the differential-algebraic functional form 

D 3 u = -2D 2 uD x u. 

Since the Lax representation for the dynamical system (II . T[) can be interpreted [1] [11] as the 
existence of a finite-dimensional invariant ideal T{u} C IC{u} realizing the corresponding finite- 
dimensional representation of the the Lie-algebraic commutator relationship p.lj) , this ideal can 
be constructed as 

(3.2) l{u} := {A Vi + Av/ 2 + z l ' 2 h € >C{u} : G K, 1 < j < 3, A G R}, 
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where v = D t u, z = D 2 u and A £ R is an arbitrary real parameter. To find finite-dimensional 
representations of the D x - and /^-differentiations, it is necessary [T] first to find the D t -invariant 
kernel ker D t C l{u] and next to check its invariance with respect to the /^-differentiation. It is 
easy to show that 

(3.3) ker D t = {/ £ IC 3 {u} : D t f = q(X)f, A £ R}, 

where the matrix q(X) := q[u, v, z; A] €End IC{u} 3 is given as 



(3.4) 



?(A) 












-A 











— Xz x 





To obtain the corresponding representation of the /^-differentiation in the space IC 3 , it suffices 
to find a matrix l(X) := l[u, v, z; A] GEnd K.{u} 3 that 

(3.5) D x f = 1(A)/ 
for / £ K,{u} 3 and the related ideal 

(3.6) K{u} := {< g,f> K3 : f £ ker A c IC 3 {u}, g £ IC 3 } 

is D x - invariant with respect to the matrfix differentiation representation 1375]) . Straightforward 
calculations using this invariance condition then yield the following matrix 



(3.7) 



/(A) 



z 

-Xtz 



\ 



A 2 v^(t 



tv) 



X 2 U\f~Z XVyfz 

—X 3 tuyfz —X 2 tv^/z 
X 4 (tuv — u 2 )— —Xv x /y r z+ 
— X 2 u x j yfz +X 3 (tv 2 — uv) —z x /2z 

entering the linear equation I3.5[) . Thus, the following proposition is proved. 

Proposition 3.1. The generalized Riemann type dynamical system U.6}) is a bi-Hamiltonian 
integrable flow possessing a non- autonomous Lax representation of the form 



\ 
) 



(3. 












-A 











— Xz x 


u x 



D x f 



—X 3 tuyfz 
X 4 (tuv — u 2 ) — 
-X 2 u x Ufz 



Xv^fz z 
-X 2 tvsfz —Xtz 
-Xv x /y/z+ X 2 ^/z(u — tv) — 



-X 3 (tv 2 



-z x /2z 



f, 



J 



where X el is an arbitrary spectral parameter and f £ C'- 00 ) (R 2 ; R 3 ). 

Remark 3.2. Simple analogs of the above differential-algebraic calculations for the case N = 2 lead 
readily to the corresponding Riemann type hydrodynamic system 

(3.9) D t u = z 2 xl D t z = 

on the functional manifold M 2 , which possesses the following matrix Lax representation: 

?- ° ), D*f=( w " X ,~ -V 
-Xz x u x J y — X(u + u x /z x —z xx /z x 

where A £ R is an arbitrary spectral parameter and / £ C^°°^ (R 2 ; R 2 ). 



(3.10) 



D t f = 



f. 



As one can readily see, these differential- algebraic results provide a direct proof of Proposition 
11.11 describing the integrability of system (|1.7[) for N — 3. The matrices (|3. 7[) are not of standard 
form since they depend explicitly on the temporal evolution parameter t £ R. Nonetheless, the 
matrices (|3. 41) and (|3.7[) satisfy for all A £ R the well-known Zakharov-Shabat type compatibility 
condition 

(3.11) D t l{X) = [ g (A), l(X)} + DJ(X) - u x l(X), 
which follows from the Lax type relationships (|3.3p and (|3.5p 

(3.12) D t f - q(X)f, D x f = l(X)f 

and the commutator condition (|3.ip . Moreover, taking into account that the dynamical system (|1.7I) 
has a compatible Poissonian pair (|2.11l) and (|2.2ip depending only on the variables (it, v, z) J £ M 3 
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and not depending on the temporal variable isK, one can certainly assume that it also possesses 
a standard autonomous Lax representation, which can possibly be found by means of a suitable 
gauge transformation of (13.121) . We plan to pursue this line of analysis in a forthcoming paper. 

4. Concluding remarks 

A new nonlinear Hamiltonian dynamical system representing Riemann type hydrodynamic equa- 
tion (jl.lj) in two and three dimensions proves to be a very interesting example of a Lax integrable 
dynamical system, as we have proved here. In particular, the integrability prerequisites of this 
dynamical system, such as compatible Poissonian structures, an infinite hierarchy of conserva- 
tion laws and related Lax representation have been constructed by means of both the symplectic 
gradient-holonomic approach [8j [TTJ [12] and innovative differential- algebraic tools devised recently 
[TJ |3] for analyzing the integrability of a special infinite hierarchy of Riemann type hydrodynamic 
systems. It is also quite clear from recent research in this area and our work in this paper that the 
dynamical system is a Lax integrable bi-Hamiltonian flow for arbitrary integers N £ N; this 
is perhaps most readily verified by means of the differential-algebraic approach, which was devised 
and successfully applied here for the cases N = 2 and 3. 

We have seen in the course of this investigation that perhaps the most important lesson that 
one can derive from this approach is the following: If an investigation of a given nonlinear Hamil- 
tonian dynamical system via the gradient-holonomic method indicates (but does not necessarily 
prove) that the system is Lax integrable, then its Lax representation, can often be shown to exist 
and then successfully derived by means of a suitably constructed invariant differential ideal I{u} 
of the ring JC{u} in accordance with the differential-algebraic approach developed here for the 
integrability analysis of the Riemann hydrodynamical system investigated above. Consequently, 
when it comes applying this lesson to the investigation of other nonlinear dynamical systems, it 
is natural to start with systems that are known to be Lax integrable and to try to identify and 
characterize those algebraic structures responsible for the existence of a related finite-dimensional 
matrix representation for the basic D x - and Dt-differentiations in a vector space K, p for some finite 

P ez+. 

It seems plausible that if one could do this for several classes of Lax integrable dynamical sys- 
tems, certain patterns in the algebraic structures may be detected that can be used to assemble a 
more extensive array of symplectic and differential-algebraic tools capable of resolving the question 
of complete integrability for many other types of nonlinear Hamiltonian dynamical systems. More- 
over, if the integrability is established in this manner, the approach should also serve as a means of 
constructing associated artifacts of the integrability such as Lax representations and hierarchies of 
mutually commuting invariants. As a particular differential- algebraic problem of interest concern- 
ing these matrix representations, one can seek to develop a scheme for the effective construction of 
functional generators of the corresponding invariant finite-dimensional ideals 1{u} C K,{u\ under 
given differential-algebraic constraints imposed on the D x - and /^-differentiations. 

We have demonstrated here that an approach combining the gradient-holonomic method with 
some recently devised differential- algebraic techniques can be a very effective and efficient way 
of investigating integrability for a particular class of infinite-dimensional Hamiltonian dynamical 
systems (generalized Riemann hydrodynamical systems) . But a closer look at the specific details of 
the approach employed here reveals, we believe, that this combination of methods can be adapted 
to perform effective integrability analysis of a much wider range of dynamical systems - a goal that 
we intend to pursue in the very near future. 
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